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^ I Based on the observation that Ko{l4:30) is lighter than its SU3 counterpart, ao(1450), we examine 

00 . the possibility that these particles, together with /o(1370), /o(1500) and /o(1710), fill a tetraquark 

recurrence of the sub-GeV 0++ nonet mixed with a glueball state. We find the picture to be consistent 
with the known data about the three /o resonances, more than the qq hypothesis. Conventional spin- 
f~^ ' orbit coupling suggests the qq, P-wave, nonet to lie around 1200 MeV. We review possible experimental 

\^ , indications of a scalar isovector resonance at 1.29 GeV, first observed by OBELIX in pp annihilation. 
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r-| ! 1 Introduction 

Meson spectroscopy has provided important clues to the understanding of the strong interactions. QCD 
^ ' is definitely well established as the basic theory, but we are still far from deriving the spectrum of mesons 
j^ ' from first principles, including the determination of which quark configurations produce reasonably de- 
fined bound states and which do not. In this area, meson phenomenology can still provide valuable input 
for more fundamental considerations. 

In this paper, we analyze the positive parity scalar mesons in the mass range 1-2 GeV. In the current 
picture, these states are supposed to arise from P-wave, qq configurations, described long ago in p. The 
well known, J^'^ = 2++, 1++, 1+" nonets, identified by a2(1320), ai(1260) and 6i(1235) respectively [^ 
do fit in the qq picture quite well. The case of 0"''+ mesons is more complicated. The well-established, 
sub-GeV scalar mesons ao(980) and /o(980) do not fit in the {qq)p-wave picture. In addition, lattice QCD 
predicts the appearence of a 0++ glueball at a mass 1500 — 1700 GeV j3], followed by 2++ and 1"' 
states at considerably higher mass. 

There are ten established scalar mesons in the 1-2 GeV range, namely: ao(1450), Kq(14:30), with 
/ = 1, 1/2 respectively, and three isosinglets: /o(1370), /o(1500) and /o(1710). It could be natural to 
identify these states with the qq nonet with L = l, 5 = 1, J = {L, S, and J denote orbital, spin and 
angular momentum of the qq system). This classification is however disfavored by two facts. 

i) Masses are too high. For S = 1, the spin-orbit coupling gives the ordering ^: m{2~^^) > m(l++) > 
to(0+"'"). Tensor forces may distort the regular spacing between multiplets, but are unlikely to push the 



scalar state on the top. In the case of charnioniuni states, for example, the scalar xo is the lowest of 
X-states, with tensor forces pushing it further down with respect to the spin-orbit coupling prediction. 

ii) Most important, i^o(1430) is lighter than ao(1450) by some 62 MeV, repeting the pattern of the 
/t(800) with respect to ao(980) and at variance with the level ordering in the other positive parity nonets^. 

Our starting point is that the reversed-octet pattern is natural for a tetraquark, [qq] [qq] nonet with 
each pair antisymmetric in flavour. Such a structure has been advocated long ago for the sub-GeV, scalar 
nonet |5, and was recently re-proposed [2] in connection with new evidence for the /o(600) and k(800) 
resonances in D non leptonic decays [5] [7], see also Ref. [S] for recent theoretical developments. 

In the present paper, we propose that the 0"*"^ decimet is made by a tetraquark recurrence of the 
light 0^^ nonet, with isosinglets mixed with the lightest glueball. 

We use a QCD inspired form of the mass matrix and of the decay amplitudes and neglect mixing with 
the lowest scalars. We find that the picture is consistent with the known information about the three /o 
resonances. We confirm the hypothesis ^1] that the /o(1500) has a dominant glueball component. As a 
byproduct, we predict the mass of the still missing, J = 1'^ , singlet in the range 1470 — 1500 MeV. 

A crucial prediction of our model is that there should be an additional nonet, the genuine qq, P-wave 
scalar, most naturally around or below the ai(1260). In this region, an additional 1 = 1 scalar state 
decaying in K^Kg was indicated time ago by the OBELIX Collaboration |3] and is not excluded by a 
more recent analysis of the rjTr resonances in p — p annihilation ^01 • There exists also some evidence for 
an additional isosinglet, e.g. the hypothetical X(lllO) T!^. It seems fair to say that the possibility of a 
nonet around or below the ai(1260) cannot be excluded at the moment. 

The problem of scalar mesons has been addressed by several authors with very different proposals, 
see |12) for a recent review. Here, we recall the investigations more closely related to our proposal. 

In the framework of the qq picture, issues (i) and (ii) above have been addressed by Schecter and 
collaborators J13j . who invoke a mixing of the higher and lower scalar nonets. The /o(1500) has been 
identified with a glueball mixed with the P-wave qq states /o(1370) and /o(1710) by Close and collab- 
orators m]. Spectrum and mixing of scalar mesons with glueballs has been studied with QCD sum 
rules since the seminal investigation in Ref. JH]- In Ref. ^j a situation similar to that of Ref. J14j is 
found, namely a glueball mixed at 20% with a qq meson^. However, QCD sum rules give also a very 
light glueball, in the range of the /o(600). Finally, the authors of Ref. ^\ argue that a single qq state 
could give rise to multiple resonances, e.g. to both ao(980) and ao(1450), due to dressing by hadronic 
interactions . 

The variety of the proposed solutions points to the need of more accurate experiments to determine 
the actual number of scalar nonets as well as their decay properties, such as will be provided for p — p 
annihilation by the new antiproton facility at GSI. 

The plan of the paper is as follows. 

We give in Sect. |21 the effective couplings used to analyse masses and decays of the mesons with 
negative parity nonets (identified by n and p), positive parity nonets (identified by oo(980), ai(1260), 
6i(1235) and a2(1320)) and the scalar decimet associated to ao(1450). We discuss the resulting mass 
spectra in Sect. |51 and decay rates in Sect.|ll where we compare our results with Refs. ^H] and ^3], and 
Sect. El We comment on the possible evidence for an additional scalar nonet in Sect.El conclusions are 
given in Sect. [71 

'^In this paper we adopt the name convention of Particle Data Book; the number used to identify each particle deviates in 
some cases from the actual value of its mass, notably: m[iCo(1430)] = (1412 ± 6) MeV, m[ao(1450)] = (1474 ± 19) MeV H]. 

^QCD sum rules have been applied to the correlation functions of pure glue operators and quark bilinears. It would be 
interesting to consider quark quadrilinears as well, to include tetraquark mesons in the picture. 



2 Effective lagrangians 

We describe nonet fields with the 3x3 matrix, 11: 
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The states associated to Xi^2 correspond to the eigenstates of the quark mass matrix: 

qq : Xi = \nn >; X2 = \ss > 

[qq] [qq] : Xi = \ [ns] [ns] >; X2 = \ [nn] [nn] >; (2) 

In the notation of Ref. ^1], n represents u and d quarks, the isospin zero combination in nn and nn is 
understood. 

We write the physical mass eigenstates according to: 

/o,iow = cos 6X1 + sin 6X2 

/o,high = - sin 9X1 + cos 9X2 (3) 

Interesting situations correspond to: ^ ~ {qq, ideal mixing); 9 ~ 90'^ ([9?] [9^]) ideal mixing); or to: 
^~-540(/o,io„~/8). 

We parametrize the nonet mass matrix according to (A = diag(0,0, (^)): 

Cm^s = m{a) + TY (A • n2) + c[Tr(n)]2 (4) 

The second term represents the effect of the strange-to-light quark mass difference. The coefficient c 
embodies the contribution of annihilation diagrams, with one quark pair transformed into any other one, 
via intermediate states made by pure glue. If S were dominant, the mass matrix would be diagonal in 
the Xi^2 basis, with c dominant it would be diagonal in the octet-singlet basis. 

We expect annihilation to give a rather small contribution in the 1-2 GeV region, therefore we neglect 
SU3 breaking in it. The mass formula Eq. @ is thus more restricted than required by a first order 
expansion in SU3 breaking. 

In the presence of a 0+^ glueball we include additional terms mixing Xi^2 with the glueball as well 
as a diagonal glueball mass, niQ = ni{a) + Anic- We neglect SU3 breaking in the annihilation diagrams 
leading from quark matter to pure glue states, as before, so that the Xi 2 to glueball mixing is described 
by a single parameter, w. In the Xi,X2, G basis, we are led to the 3x3 matrix: 

/ 2c V2c V^w \ 
M = m{a) + V2c c + 26 w \ (5) 

\ \/2u; w AniG / 

We shall fix the combination S = 3c -|- Anic from the sum of isosinglet masses and take 5 from the 
Kq — oq mass difference. Setting, in addition, Amc = 5^(1 — a), we arrive at the 2-parameter formula: 

M = m{a) + I ^^Q,S iaS + 26 w \ (6) 

V2w w (1 — a)S 



with: S = X^j '72-(/o ) — 25 — ?>m{a), a and w independent parameters. The same formula with a = 
was proposed in Ref. [Tsj . 

Next we consider decay rates. 

We restrict to decays into pairs of pseudoscalar, strictly octet, states. Amplitudes for the decay in 
isosinglet particles, i.e. r] or ry', involve additional, undetermined parameters. In the exact SU^ limit, the 
effective lagrangian contains three terms: 

Cdec^y = ^{[Tr(n • P2) - X Tr(n)Tr(p2)] + 7 G TV(p2)} (7) 

where P is the pseudoscalar meson matrix, the analog of (^, G the glueball field, A, x and 7 unknown 
parameters, the latter describing the relative coupling of the glueball to two meson states. 

As a further simplification, we require that the quarks originally present in the decaying meson appear 
in the final mesons. It is an approximation related to the OZI rule |19j which forbids, e.g. (/)(1020) — > tttt, 
and is related to the smallness of annihilation processes. The rule restricts the value of x by requiring 
that either X2 (for qq) or Xi (for [(/(/][g^]) do not decay to vrvr. This gives: 

qq : X = (8) 

or 

[qq][qq]: x = -1/2 (9) 

For singlet decays, we have to substitute the combinations corresponding to the eigenstates of ©, 
namely: 

f^'^=cfXi + cfX2 + cfG- (i = 1,2,3) (10) 

If we call Qh the coefficient in (0 of a given decay channel, the decay rate is: 

- - 'P' \F{v')\'9l (11) 



" 87r3M2 

M is the mass of the decaying particle, p the decay momentum and -F(p^) a form factor suppressing 
the emission of fast particles by the initial meson, due to its extended structure. In Ref. jjl], the 

parametrization is used: 

|F(p2)|2 = e-pV8/3^ (12) 

with (3 = 400 MeV. We shall consider also the extreme cases of no form factor (/? = +00) or /? = 240 MeV, 
to exemplify a more extended structure than qq. 

3 The mass spectrum 

For each nonet, the parameter 5 is obtained from K — a mass splitting and c can be fixed by the value of 
the sum of isosinglet masses. Thus, we can predict the difference of the isosinglet masses, Amj, as well 
as the mixing angle. 

To test the accuracy of Eq. Q , we have reported in Tables ^ and [21 the predicted values of Amj 
compared to the experimental ones, given in parentheses, for the well identified pseudoscalar, vector, 
tensor and axial vector nonets and for the two scalar multiplets. The 5 parameters and the Xi — X2 
mixing angles are also reported. 

Masses have been taken from ,2. except for the lowest lying 0"*"^ nonet, for which we have taken the 
values quoted in Ref. 0. The values quoted by the BES Collaboration [7] are consistent, within the large 
experimental errors. 
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Table 1: Multiplets below 1 GeV. Above, for linear mass formulae: i) predicted values of the isosinglets 
mass difference, Amj, compared to the experimental values (in parentheses); ii) SU3 breaking parameter, 
6; iii) isosinglet mixing angles. Below: same for quadratic mass formulae. 



Our mass formula can be applied with squared masses. Quadratic mass formulae work somewhat 
better, as seen from the Tables, in particular for the pseudoscalar mesons. For simplicity of exposition, 
we shall restrict to linear mass formulae in what follows. 

The experimental spectra are illustrated in Fig. ^tolSl The difference of the pattern exhibited by the 
scalar multiplets sticks out conspicuously and is reflected in the negative sign of 5 in Tables H and El 

A few comments are in order. 

• With the well-known exception of the pseudoscalar mesons, the predicted values of Amj are in 
general agreement with observation, and mixing angles are all close to maximal, Eq. Q, thus 
supporting the Ansatz made for the mass formula; 

• Mixing angles reproduce well the affinities for non-strange (vrvr, /jvr, etc.) or strange {KK, K*K, 
etc.) decay channels of the qq nonets and of the lowest scalar nonet, the latter under the [(7(7] [gg] 
hypothesis |^; 

• The contribution of gluon annihilation to the mass matrix is rather small as expected. As an 
example, we find for the the tensor mesons: 
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Figure 1: Mass spectra of the lowest lying meson multiplets, pseudoscalar, vector and the scalar mesons below 1 GeV. In 
blue the 1=1 states, pink and orange / = 1/2 and / = respectively. Masses and particle symbols are taken from Particle 
Data Group, Ref. [5], except for /o(600) and fi:(800) where we use the values given by E791 0. 
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Table 2: Same as in Table ^ for the multiplets in 1-2 GeV. For the 1"*" the experimental Aruf is not 
available since /ig has not yet been observed. 

• The two axial vector nonets differ in charge conjugation but K-states are mixed by SU^ breaking. 
The mass matrix contains seven parameters which can be determined from the known masses, 
yielding a prediction for the yet unknown h'l , with a slight difference between linear and quadratic 
mass formula: 



m{h[) = 1474 MeV (quadratic); m{h[) = 1498 MeV (linear) 



(14) 



To fit the masses, we need a rather small non diagonal element of the Hamiltonian between Kaon 
states: 

< K++\H\K+^ >= 53 MeV (15) 

The results in (|13() and H15() indicate rather small configuration mixing and annihilation contributions 
in the meson states above one GeV. 
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Figure 2: Mass spectra of the axial mesons multiplets, same conventions as in Fig. Strange states are mixed by SU3 
breaking and are reported separately. The green-dashed line indicates the predicted position of h'l , yet to be observed state. 

For the O^"*" decimet, we solve numerically the eigenvalue equation corresponding to Eq. ©. We use 
masses and errors given in j2] to compute the x^ of each determination {x^ = J2i(.''^th — ''^exp)'^ / ^i^l^p) 
except for /o(1370) whose mass is rather poorly known. In this case we take, for definiteness: m = 
1370 ± 60 MeV. 

We show in Fig. 0]the distribution of x^ iii parameter space^. The red region gives a better x'^^ but 
acceptable masses are found in all the green region. In the left corner (w ~ 0, q small) the glueball 
corresponds to the highest mass and it moves to the intermediate mass at the right corner (tD ~ 0, a 



^w can be made to be positive by an appropriate definition of the sign of the G field 
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Figure 3: Mass spectra of the tensor mesons and of the higher scalar mesons. Same conventions of the previous figures. 
Note the inversion of the ordering of J = 1 vs. / — 1/2 states for scalar mesons, with respect to the well established 
i<l<j)s/p~wave multiplets. In addition, there is an extra 0"*"^, isoscalar, state, indicated by the green-dashed line. 



large). Of the other two states, the highest one is mostly Xi and the lowest X2- For qq states this 
corresponds to the ss state being lighter than the nn and it is, of course, a consequence of the unnatural 
sign of the effective strange-to-light quark mass difference, implied by the Kq—uq splitting. For tetraquark 
states, on the contrary, this is the natural level ordering, Eq. 0. 
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Figure 4: Distribution in parameter space of the x^ computed from the eigenvalues of the mass matrix Q. The arrow in 
the right panel points to the region preferred by particle decays in the [gg][gg] hypothesis. The bars in the panel give the 
probabilities of quark/glueball composition of the eigenvectors. The left panel corresponds to the qq hypothesis. 



4 Decay rates of singlet mesons and best solution 

For singlet meson decays, we have explored the parameter space by optimizing the x^ with respect to 
the ratios of branching fractions reported in the first two columns of Table El 

In the [g(7][gg] hypothesis, a reasonable picture is obtained for the point indicated in Fig. [fright panel, 
with: a = 0.81, u;/S = 0.16 and 7 = -0.26. 

For reference, we give explicitly the mass matrix (|16|) . in the basis Xi, X2, G: 



1626 108 64 
M=\ 108 1426 45 
64 45 1528 



(16) 



and the corresponding eigenvectors and eigenvalues: 
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(17) 



The value mc = 1528 is consistent with lattice determinations P). 

The numerical values of the branching fraction ratios computed in correspondence to Eq. (4qmix) are 
given in Table IHl Values generally agree with the trend of the experimental figures. Larger form factor 
damping is preferred. Values in the third column reproduce well the affinity of /o(1370) to non strange 
channels (indicated independently by the large branching ratio in four pions: i?47r[/o(1370)] ~ 80% |2]). 

The same procedure with the qq hypothesis gives as best solution the point indicated by the left panel 
in Fig. ^ corresponding to a = 0.23, w/Ti = 0.14. The mass matrix and corresponding eigenvalues and 
eigenvectors for this case are: 



1519 32 56 
M=\ 32 1373 39 
56 39 1688 



(18) 
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Branching fraction ratios for the qq assumption are also reported in Table IHI For the first two columns 
we find reasonable agreement with data. However, within our model of the mass matrix, the lowest state 
is essentially ss, see Fig. and the /o(1370) is predicted to have a strong affinity for strange decay 
channels, in contrast with data. 
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Table 3: Ratios of branching fractions for vrvr and KK decays of the 0^^ isosinglets. Above: [(7(7] [^^] 
hypothesis with three different assumptions on the form factor. Middle: qq hypothesis (no form factor 
gives the best agreement). Below: experimental data from Particle Listing 2005 [2]: a |20j . b |21j . c |22j . 
d I2F 



Before closing this Section, we compare our results with the previous proposals in [T3| and J14l . 

In ^3]) the anomalous properties of the O"*"*" multiplet above 1 GeV are ascribed to a mixing of the 
P-wave qq states with the the sub-GeV scalar mesons, assumed to be of a different kind (i.e. hadronic 
molecules or tetraquarks) . A suitable mixing of the two multiplets could push the qq scalars above what 
expected in the spin-orbit coupling scheme and in addition make the 1 = 1 state heavier than the / = 1/2 
one. 



The mixing needs to be very fine-tuned and, in addition, it turns out to be quite large. Quantitatively, 
for linear mass-formulae, a non diagonal matrix element of the Hamiltonian of about 250 MeV is required 
in Ref. ^S] to push the oo(1450) above the Kq{1430). The result is rather unnatural given the very differ- 
ent configurations of the states involved {qq vs. qqqq) and that it involves quark pair annihilation. In the 
more favourable case of the strange mesons belonging to the 1^^ and 1^ octets (both qq configuration, 
no annihilation required) we have found a quite smaller result, Eq. (|15() . On our side, the non diagonal 
matrix elements in Eqs. (|16[ I18 |) . are similar to (|15|) . 

Comparison with ^^ requires some qualification. Without mixing to the sub-GeV multiplet, the 
inverted pattern of the 1 = 1 and / = 1/2 masses could still be explained, in the q — q model, as due 
to an overcompensation of the strange to non strange quark mass difference by spin effects (which also 
depend upon quark masses). Within our Ansatz for the mass matrix, however, the wrong sign of the SU3 
breaking in the octet propagates to the singlets and makes it so that the scalar made by strange quarks 
is lighter than the one made by non strange quarks. For this reason, the gg+glueball assumption fails in 
reproducing the affinity of /o(1370) for non strange channels. Table IHI The tetraquark+glueball model 
fares better in that the association of the heavier meson with strange quarks is preserved in the reversed 
octet. 

With a more general mass matrix, the correlation between octet and singlets would be lost and the 
observed affinities to non strange or strange channels could be reproduced. The difficulty would remain, 
however, to explain the opposite sign of SU3 breaking in the O"*""*" octet with respect to the other P-wave 
octets. 

5 The decay rates: ao(1450) -^ KK and 7^0(1430) -^ Kit 

In the exact SU3 limit, these partial decays are given in terms of one single amplitude. Comparing with 
data, we can test wether the two particles belong indeed to the same SU3 octet. Ref. [2j quotes: 

r[iv:o(1430)] = 294 ± 23 MeV 

B[Ko{U30) -^ Kit] = 0.93 ± 0.10 (20) 

whereby: 

r[i^o(1430) -^ Ktt] = 253 ± 51 MeV (21) 

Using (|TT|) . in the SU3 limit we predict: 

9 A/f^ 

r(a+ ^ KK)sus = o — 7#r(i^o+ ^ Ktt) = 150 ± 28 MeV (22) 

As for data on ao(1450) decays, Ref. |2j quotes three observed decay channels: ryvr, 77' vr and ujtttt, and 
gives total width and ratios of the partial widths to the rjir one: 

r(ao(1450)) = 265 ± 13 MeV; 

r(r?'vr) T(KK) T(ujit-k) . , , , . x 

^V , = 0.35 ± 0.16; -^7 — r^ = 0.88 ± 0.23; -^7 — -^ = preferred value not quoted (23) 
r(777r) r(r/7r) '^[W) 

Neglecting the lotttt contribution, we find: 

r(a+ -^ KK) = 105 ± 46 MeV (24) 

9 



In ^^ the discrepancy of the central values was taken to support unequal mixing angles of ao(1450) 
and /fo(1430) with their lighter partners. However, this would be an SU3 breaking effect of the same 
order of similar effects neglected in decay rates, so the argument may not be so compelling. In fact, given 
the large experimental errors and the neglect of SU3 breaking, the assignment to the same multiplet can 
be consistent. 

A reduction of the experimental errors and, most important, a determination of the partial width of 
Oo(1450) — > WTTvr would greatly help to clarify the situation. 

6 Another scalar nonet? 

The crucial feature of our model is that it requires an additional nonet, the genuine qq, P-wave scalar, 
around or below the oi(1260). Experimental searches for SU3 scalar nonets have been focussed on 
isovector resonances decaying into K^Kg or rjir. 

The very accurate analysis oi p — p annihilation at rest done by the OBELIX Collaboration 9J showed 
a clear peak in the mass distribution of K^Kg around 1.3 GeV, where the 02(1320) is expected. However, 
a comparison of annihilation events produced in targets with different density showed that the peak shifts 
from 1.32 GeV at the lowest density target (where most of the reaction proceeds from P-wave), to 1.29 
GeV in the liquid Hydrogen target (where annihilation from S-wave dominates). The best fit for the 1.29 
GeV resonance is obtained for J^*-^ = O"*""*", indicating an additional scalar isovector ao(1290), produced 
dominantly by S-wave annihilation. 

The most recent, high-statistics, analysis presented in JHl used p — p annihilation in flight, where 
it is not possible to separate the different partial waves, and the t]tt channel. Also in this case, a peak 
around 1.3 GeV was seen. The best fit was indeed obtained for J = 0"'"'', with mass consistent with 
ao(1290) but a larger width than indicated by OBELIX. A lower quality fit is obtained for the 02(1320), 
which was preferred "because of the weight of numerous final states, including ryvr , in which O2(1320) 
has been observed, and not ao(1290)". In view of the OBELIX data, the most likely situation would be 
in fact that the two resonances coexist in the final states. This hypothesis was also tried in jJOl but with 
inconclusive results. 

It seems to us fair to say that the possibility of a nonet around or below the ai(1260) cannot be 
excluded at the moment and that the issue of the oo(1290), a good candidate for the qq P-wave scalar, 
deserves a closer scrutiny. 

7 Conclusions 

Scalar mesons present many puzzling aspects, which have given rise to several different speculations. 

We have presented a novel proposal in which the presently observed scalars, both below 1 GeV and 
in the 1-2 GeV range, are tetraquark states. This explain naturally the mass ordering of / = 1/2 and 
1=1 states exhibited by both multiplets and is consistent with known facts about their decays. The 
/o(1500) is naturally included in the picture as a light glueball, mixed with the tetraquark singlet states 
/o(1370) and /o(1710). 

A direct consequence of our proposal is that there should exist another scalar nonet, probably around 
1.2 GeV, the genuine qq P-wave state. Evidence of such a complex may have been seen by OBELIX in pp 
annihilation at rest and may be indicated, or anyway not excluded, by a recent high statistics experiment 
on pp annihilation at higher energy. 
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A thorough exploration of the scalar spectrum in the 1-2 GeV region and a more precise determination 
of decay modes and branching ratios is clearly needed for the resolution of the nature of scalar mesons, 
an issue of the utmost importance for understanding the way QCD works for bound states. 

Acknowledgments. We thank F. Close for interesting exchange on the properties of the /o(1500) and 
U. Gastaldi, A. Masoni and C. Voena for useful information on the experimental data. 
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